Landau Model for Commensurate-Commensurate Phase Transitions in Uniaxial 

Improper Ferroelectric Crystals 



C/3 



o 



> 



M. Latkovic and A. Bjelis 
Department of Theoretical Physics, Faculty of Science 
University of Zagreb, Bijenicka 32, 10000 Zagreb, Croatia 

V. Dananic 



o 
o 

o . 

, Department of Physics, Faculty of Chemical Engineering and Technology 



University of Zagreb, Marulicev trg 19, 10000 Zagreb, Croatia 



^ We propose the Landau model for lock-in phase transitions in uniaxially modulated improper 

'^'h [ ferroelectric incommensurate-commensurate systems of class I. It includes Umklapp terms of third 

and fourth order and secondary order parameter representing the local polarization. The corre- 
sponding phase diagram has the structure of harmless staircase, with the allowed wave numbers 
obeying the Farey tree algorithm. Among the stable commensurate phases only those with periods 
equal to odd number of lattice constants have finite macroscopic polarizations. These results are in 
^ ' excellent agreement with experimental findings in some A2BX4 compounds. 
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Many uniaxial systems undergo a sequences of incommensurate-commensurate (IC-C) or commensurate- 
commensurate (C-C) phase transitions It is generally accepted that the main cause of such behavior comes 
"j^ from the mode softening induced by competing interactions between neighboring sites in the crystal lattice . The 
mode softening occurring at an arbitrary point in the Brillouin zone characterizes the so called class I of IC-C systems, 
I \ in contrast to the examples of the class II for which the wave vectors of soft modes are at (or close to) the center or 
■ O ■ the border of the Brillouin zone |^ . 

The IC-C and C-C phase transitions are mostly of the first order, and are accompanied with hysteresis and memory 
effects. Also, depending on the wave number of modulation, some commensurate phases have a uniform (e.g. ferro- 
electric) component. Recent experiments show that the modulation of commensurate phases is usually domain like, 
suggesting that multisoliton configurations that are usually present in the incommensurate states, are to some extent 
frozen in lock-in phases as more or less dense soliton lattices [|]-|7|. 
\Q ■ Theoretical considerations of the class I of IC-C systems mostly start either from discrete microscopic models 
\0 ' of competing interactions or from phenomenological expansions of the free energy functional. In the former cases 
i the sequence of phase transitions is usually characterized by complete or incomplete devil's staircases for the wave 
">! I number of ordering . The crucial importance of discreteness in these models comes from the assumption that 

■ the couplings between neighbouring atoms or molecules are very strong. However, as argued below, this assumption 
_ usually does not suit the microscopic properties of real materials. 

■ The latter approach is based on the Landau theory of phase transitions, and is generally appropriate for weakly 
Ci ' coupled systems. The justification for continuous Landau models comes from many experimental indications, e.g. 
^ from the neutron scattering data |pT|-|l^ , showing a well-defined dispersion curves for collective modes with distinct 

' soft-mode minima. The basic Landau model, commonly applied to various systems of the class I, includes one Umklapp 
^ i term and the Lifshitz invariant. Within the mean field approximation it leads to the sine-Gordon problem |l^-|l^] 
Q I which has one isolated phase transition of the second order to the unique commensurate phase. However, in such 
O ' systems one usually encounters more commensurate phases, sometimes with the order of commensurability much 
^ ■ higher than three or four. 
• The phases with high order commensurabilities can be explained only by some extension of the sine-Gordon model. 

/\ One widely explored way is to introduce an additional set of high order Umklapp terms p^-|2C[| , more precisely to take 
^ into account as many terms as there are stable commensurate states in the experimental phase diagram. This approach 
■ " " is however not free from serious disadvantages. At first, one introduces terms with powers of the order parameter 
amplitude which are twice the order of commensurability, and as such are beyond the standard (minimal) Landau 
scheme that includes only those higher (usually fourth) order terms which necessarily guarantee the boundness of the 
free energy density from below. Further on and even more important, the microscopic analyses, which are usually 
avoided in such approaches, show that such terms are as a rule very small, and as such physically irrelevant in weakly 
coupled systems. 

The model which is, in contrast to those of Refs. |p^-po|, formulated within the standard Landau scheme, and still 
explains the stabilization of high order commensurate phases, is proposed recently by two of us [Elj,E2|. In this model 
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we take into considerations two Umklapp terms of the lowest possible orders in the free energy expansion, i.e. those of 
the third and fourth order, and so remain strictly within the Landau theory |2^. The similar starting point, but with 



Umklapp terms of orders higher than four, was proposed earlier in Ref. |g4|. In Refs. |21|J2^ we show that, due to 
the competition between these two Umklapp terms, the periods of ordering follow harmless staircase with the values 
obeying a Farey tree algorithm, while the corresponding commensurate configurations have the soliton lattice-like 
forms. Also, the hysteresis and memory effects which are usually observed in particular materials, are within the 
present model interpreted as an intrinsic property, i.e. as consequences of the free energy barriers that appear due 
to the nonintegrability of the Landau functional. These barriers prevent smooth transitions between neighbouring 
thermodynamically stable configurations |p5|- p7| , ^ , p2[ . 

The Landau model from Refs. [ pTp^ contains only one order parameter, and as such does not possess all ingredients 
needed for the explanation of phase transitions in uniaxial improper ferroelectric materials. For such materials it is 
necessary to include a secondary order parameter, that responsible for the net polarization of particular modulation 
structures. Although the necessity for such extension was put forward already by lizumi et al. the subsequent 
studies of phase transitions in improper ferroelectrics were developed mostly within the simple sine-Gordon model 
with one Umklapp term ||2^,^ . These analyses were concentrated on one isolated IC-C phase transition, allowing at 
most for its first order nature. 

Our aim is to show that the model of lizumi et al. extended with the additional Umklapp term which provides 
the mechanism for sequences of C-C phase transitions [ ^p2| , enables the understanding of the complex phase diagrams 
in A2BX4 compounds. In particular we explain several features of C-C phase transitions, like harmless staircase 
behaviour of the wave number of ordering, the first order nature of phase transitions, and the polarization and 
modulation properties of commensurate phases. 

We start from the assumption that the quadratic soft mode contribution to the Landau expansion has minima at 
wave numbers {+qc, —Qc), where < qc < (73, with q^ = 27r/4 and (73 = 27r/3 (the unit length is taken equal to the 
lattice constant). The distances of qc from 173 and (74 are denoted by S3 and S4 respectively, with 53 + 84 = 7r/6, which 
gives the first independent control parameter, 64. The order parameter is complex, pe*"^. Limiting the further analysis 
to the temperature range well below the critical temperature for the transition from the normal (disordered) to the 
incommensurate phase, we also make the usual approximation of space independent amplitude p, and keep only the 
phase-dependent part of the free energy density. 

To allow for a finite polarization along some direction, we introduce a secondary order parameter u, and assume 
that it is coupled to the third-order Umklapp term, in agreement with the symmetry requirements for A2BX4 crystal 
lattice In order to stabilize the free energy with respect to changes in u, we also introduce the terms proportional 
to and u'^. The complete free energy density reads 

fi(t),u,x) = + + ^\u^ + Bu cos[3(j) + - S4)x] + C COS (4(1) - 4S4X), (1) 

where primes denote spatial derivation. The coefficients A, B and C in Eq. (^ are rescaled in order to simplify 
equations. Thus f{(f),u,x) is the original free energy density divided by p^. Note also that after this rescaling the 
coefficients B and C in Eq. (|l|) are linear and quadratic in the order parameter amplitude p respectively. Altogether, 
there are three control parameters, namely 84, B and C , while the parameter A just defines a scale for the polarization 
u. Among the elastic terms in Eq. (|l]), the first one ((/)'-dependent) favors the incommensurate sinusoidal ordering 
with the wave number q^.. The last two Umklapp terms cause the harmless staircase behaviour of the wave number 
of ordering [pi|,p2|. The above Landau expansion follows, after spatial continuation, from microscopic models like e.g. 



those of Refs. [ p0 | j31[ , that take into account local interactions presumably responsible for the mode softening (and for 
the coupling to the secondary order parameter), and start from the discrete presentation along the uniaxial direction. 
The mean-field approximation for free energy (|l|) leads to the Euler-Lagrange (EL) equations: 

(/)"+3Busin [30-K3(- ~ 54)x] + 4Csin (40 - 4(54^;) 0, 
6 

u"-Xu- Bcos[i(t> + i{- - 54)x\^Q. (2) 
6 

We are interested only in those solutions of EL equations which participate in the thermodynamic phase diagram of 
the model (^. These solutions have, for given fixed values of the control parameters, the lowest values of the averaged 
free energy 

{F) = ^ fdxf[^ix),u{x),x]. (3) 
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L is the macroscopic length of the system. 

Before calculating such solutions, let us make few remarks on EL equations (H). They can be considered as 
Lagrange equations for an equivalent classical mechanical problem. Note however that general ("classical mechanical") 
solutions of these equations are not bounded, since A is positive. The only bounded solutions are periodic. From 
one side, obviously only such solutions may participate in the thermodynamic phase diagram. From the other side, 
the unboundedness of all other solutions makes the analysis of the present model more complicated than that of the 
model without the secondary order parameter u | |2l| , p^ , for which the mechanical phase portrait is bounded, although 
chaotic 

It is clear from above remarks that periodic solutions are orbitally unstable. Therefore they cannot be calculated 
by a direct numerical integration of EL equations (^. Before embarking into another suitable numerical method, it is 
useful to establish, by extending straightforwardly the previous treatment , necessary analytic conditions for these 
periodic solutions. We start by looking for allowed periods, after taking into account that periodic solutions have to 
obey the relations 

(l){x + P) = <i){x) + (j)p, 

u{x + P) = u{x), (4) 

where P is the period and (j>p is the phase increment per period. From EL equations (H) it follows that allowed periods 
and phase increments are 

P = 4k + 3l, (l,p^S4P~l^, (5) 

where k and I are integers. The corresponding values of the total wave number (measured from the origin of Brillouin 
zone) are then 

The values of q given by Eq. (^ form a Farey tree structure, shown in Fig. |l| for wave numbers between (7=1/3 
(A: = 0, / = 1 and P = 3) and q = 1/4 {k = 1, I = and P = 4). The periodic solutions with q = 1/3 and q = 1/4 
are the basic commensurate configurations, favoured by the Umklapp terms of third and fourth order respectively. 
A higher order commensurate configurations are represented by periodic solutions with wave numbers q positioned 
between q = 1/3 and q = 1/4 in the Farey tree. The values of integers k and I are positive for all these configurations. 
We emphasize that the Farey tree structure for q is not the characteristic of usual continuous Landau models for 
IC-C transitions Jlsj-pTt. Here it is imposed through the competition between two Umklapp terms in the first of EL 
equations 

For numerical purposes it appears convenient to eliminate the explicit x— dependence from one of Umklapp terms 
in the EL equations (^, by passing from the variable (t){x) to 

ij{x) = cl>{x) + ( J - 6i)x. (7) 



The EL equations (0) now read 
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ip" + 3Bu sin (3V') + 4C sin {iip a;) = 

3 



u" - Aw - Bcos(3V) = 0. (8) 
The corresponding free energy acquires a (54-dependent term in the form of Lifshitz invariant, 

F= [ dxi-ltP' -{-^ Si)]^ + +-u''^ + -\u^ + Bu cos {3ij) + Ccos {4ij^—x)}, (9) 
J 2 6 2 2 3 

which considerably simplifies the calculation of the (54-dependence of the averaged free energy (||) for particular periodic 
solutions of the EL equations (|8|). Note that these equations do not contain the parameter 64, so that it is sufficient 
to find solutions in variable i^ix) for some parameter values of B and C. 

Once the analytic conditions (^) are established, the orbitally unstable periodic solutions of Eq. (||) can be system- 
atically calculated by treating Eqs. (j^) as a boundary value problem, and using an appropriate numerical algorithm 
suitable for nonlinear differential equations, like the finite difference method. The boundary conditions have to be 
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specified on the left and right end points of the integration, i. e. for x = Xq and x = xq + P- It is convenient to put 
X = xq at one of the inflection points of ip{x) . To complete the boundary conditions we still have to find out values 
of ^(xo) and one of the values u{xo) or u'{xo), as will be discussed in detail elsewhere Tab. | contains choices 
of Xq, and u{xq) or u'{xo) that specify completely the appropriate boundary conditions. As previously pl| , p2| , 

we find two independent periodic solutions for given values of k and I (i.e. for a given wave number q). They differ 
by symmetry and have distinct average free energies. These solutions are denoted as type A (antisymmetric in ip, 
symmetric in u) and type B (without particular symmetry). These solutions are uniquely determined once the param- 
eters {k,l) [i.e. {P,(f)p)] are chosen. Note that they do not depend continuously on the parameter S4. The additional 
minimization of the free energy (^ provides only the i54-dependence of ranges of stability for each pair (fc, I). 

Three low order periodic solutions of Eq. (||), with wave numbers q = 1/3, 1/4 and 2/7, are shown in Fig. ^ For 
each of these values the phase of the primary order parameter, ^{x), shows qualitatively different beavior. A linear 
^-dependence of tp (i.e. a simple sinusoidal modulation of the primary order parameter) is realized for q = 1/4. For 
q = 1/3 ip has an additional sinusoidal variation, while for q = 2/7 one encounters the periodic alternation of short 
domains with g = 1/4 and q = 1/3 modulations, or, in other words, the formation of a dense soliton lattice. The 
tendency towards more and more dilute soliton lattices strengthens as one goes down along the Farey treee from 
Fig. 1^. The secondary order parameter u{x), i.e. the local polarization, either changes periodically in space with 
alternate positive and negative values forming an antiferroelectric lattice (type B solution from Fig. [^a, and type 
A and B solutions from Fig. and ^c), or has everywhere a ferroelectric space dependence (type A solution from 
Fig.|a). 

From the macroscopic point of view we are interested in spatially averaged polarizations, given by 



1 



Xo + P 



(u) = — / u{x)dx (10) 

Jxo 

for periodic solutions. The values of (u) for solutions with various values of the wave number q are listed in Tab. ||. 
Note that the averaged polarization (u) vanishes for all solutions with even values, and all type B solutions with odd 
values, of the period P. 

The above boundary value method enables the calculation of thermodynamically stable solutions for rather high 
values of the parameters k and I in Eq. (|^). Here we limit the analysis of the phase diagram by keeping up to the 
tenth row in the Farey tree from Fig. |l|. Also, in order to facilitate the further discussion, we fix the value of the 
parameter B, and calculate the phase diagram in the reduced parametric space (C, ^4), as shown in Fig. |3[ Other 
choices of the cross sections in the parameter space [ p2[ lead to the qualitatively same conclusions. 

As is seen in Fig. only the lowest order commensurate phases, namely 1/3, 1/4, 2/7, and 3/10 are present 
for large values of C. By decreasing C, i.e. by moving towards the sine-Gordon limit (C —>■ 0), more and more 
higher order commensurate phases start to participate in the phase diagram. E.g., for C ~ 0.05 new phases with the 
commensurabilities 3/11, 4/13, 5/16 are present. 

The phases which have nonzero average polarizations ( p^ ) are denoted by the letter F in Fig. ^. As numerical 
calculations show, all phases with odd periods, e.g. those with wave numbers 1/3, 2/7 and 3/11, turn out to be 
ferroelectric. In other words, configurations with odd periods that have lowest averaged free energies, all belong to 
the type A solutions from Tab. H Note also that all lines in Fig. || represent phase transitions of the first order. 

The phase diagram from FigTH is in qualitative agreement with experimental phase diagrams for some members 
of the A2BX4 family. As an example we take rubidium tetrabromozincate, Rb2ZnBr4, for which there is a variety of 
data on the temperature and pressure variation of the modulation wave number pll-pJ3^,p4l-p6| . Our phase diagram 



is in qualitative agreement with high pressure measurements collected in Fig. 9 of Ref. |34| and in Fig. 8 of Ref. [ p6[ . 
In particular, the experimentally observed high order commensurate phases with q equal to 2/7, 3/10, 3/11, 4/15, 
5/17 and 7/24 are all present in Fig. ^. Note also that our expression (I) (and Fig. 0), which is obtained as an 
inherent property of the competition between two Umklapp terms |3^ , represents the theoretic explanation of the 
phenomenological hint on the Farey tree structure of experimentally observed commensurabilities p^-^ . Regarding 
the macroscopic polarization associated with various commensurate phases with odd periods, up to now the measure- 



ments are reported only for that with the lowest commensurability (1/3) 134 3^]. It shows a finite ferroelectric order, 
in agreement with our results. The experimental investigations of phases with higher commensurabilities, for which 
our results suggest presence or absence of a finite ferroelectric component, are highly desirable. 

The next example from A2BX4 family of is ammonium tetrachlorozincate, (NH4)2ZnCl4. Its pressure-temperature 
phase diagram [39| ] is similar to that of Rb2ZnBr4. In particular it contains a phase sequence with commensurabilities 
1/3, 2/7, and 1/4, with some ambiguities about region of the stability of phase 2/7 (see Ref. and references 
therein). The phase 1/3 is ferroelectric, the mixture of phases 2/7 and 1/4 shows a weak ferroelectricity, and some 
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antiferroelectric order is observed in the nonpolar phase 1/4. This is in agreement with our results in Fig. ||. Indeed, 
the local polarization u{x) for type A solution is constant for phase 1/3, it is nonpolar and has an antiferroelectric 
modulation for phase 1/4, while for the phase 2/7 we expect a weak ferroelectricity since the solution u{x) of type A, 
although almost antiferromagnetically modulated, still has a nonzero averaged value (u) . 

The last example, ammonium hydrogen selenate |41|, NH4HSe04 (together with its deuterated version ND4DSe04), 
does not belong to the A2BX4 family. The sequence of phase transitions shown in Fig. 8 of Rcf. namely 1/3 (F), 
3/10, 2/7 (F) and 1/4, can be also reproduced by going along a particular path in Fig. ^. Let us also note that the 
present model is not directly applicable to betaine-calciumchloride-dihydrate (BCCD), a system with a well-known 
rich sequence of IC-C and C-C phase transitions 43|. In this material, usually considered as a member of class 
II, one very probably encounters a competition of the first (ferroelectric, q = 0) and some higher order (probably 
q — 1/4) commensurabilities. 

Let us finally shortly focus on the form of the commensurate modulations. As already mentioned, multisoliton 
structures are sometimes frozen in commensurate phases well below the IC-C transitions Our numerical analysis 

shows that the stable solutions of higher orders (like that with q = 2/7 shown in Fig. He) have the properties of dense, 
and not dilute, soliton lattices. It also suggests that particular solutions keep this form well below lock-in phase 
transitions, and that the polarization closely follows such behaviour of the dominant order parameter by forming the 
modulated ferroelectric or antiferroelectric patterns. 

In conclusion, by extending the basic model [ pT]j2^ with the secondary order parameter representing the ferro- 
electric polarization, we obtain the phase diagram for the series of commensurate phases which is in qualitative and 
quantitative agreement with ex; 



:perimental phase diagrams for some A2BX4 compounds |34|- |3q , |39| , [40|| , as well as for 
ammonium hydrogen selenate [4l|. The main results of the present analysis are as follows. The wave number of 
ordering is given by the harmless staircase obeying the Farey tree algorithm. The commensurate phases are mainly 
characterized by dense soliton lattice modulations. The phase transitions between successive lock-in phases are of the 
first order. Finally, we establish the selection rule, by which only those among the commensurate phases which have 
odd periods are accompanied by finite macroscopic electric polarizations. 

The work is supported by the Ministry of Science and Technology of the Republic of Croatia through the project 
no. 119201. 
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TABLE I. Boundary conditions for periodic solutions of EL equations (0). 
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TABLE H. Polarization properties of commensurate phases with wave number q (Eq. Letter F indicates that particular 
solution have nonzero average polarization (Eq. ([l^)). 
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FIG. 1. Farey tree for wave numbers q defined by Eq. 
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FIG. 2. Solutions of EL equations (g) (phase of the primary order parameter) and u{x) (local polarization): k — 0, 

1^1, q = 1/3 (a), fc = 1, Z = 0, (? = 1/4 (b) and k = 1, I = 1, q = 2/7 (c). Other parameters are B = 0.1, C = 0.1, A = 1. 
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FIG. 3. Phase diagram in the (0, 54) plane for B = 0.3 and A = 1. The numbers in the figure represents wave numbers of 
ordering q of commensurate phases. Phases with nonzero average polarizations (w) are marked by the letter F. 
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